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Abstract 



This paper presents a new non-associative algebra which is used to (i) show 
how the spectator (or Gross) two-body equations and electromagnetic currents 
can be formally derived from the Bethe-Salpeter equation and currents if both 
are treated to all orders, (ii) obtain explicit expressions for the Gross two-body 
electromagnetic currents valid to any order, and (iii) prove that the currents so 
derived are exactly gauge invariant when truncated consistently to any finite 
order. In addition to presenting these new results, this work complements and 
extends previous treatments based largely on the analysis of sums of Feynman 
diagrams. 
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I. INTRODUCTION 



The two-body spectator (or Gross) equations were first introduced in 1969 and have 
been developed in a number of subsequent papers [|]]. The treatment of electromagnetic 
interactions in this context has also been studied ||- ||. However, all of these previous 
treatments have been largely based on the analysis of Feynman diagrams, and the equations 
have been largely derived from this diagrammatic analysis. In this paper we present an 
algebraic derivation of the equations which is complementary to previous diagrammatic 
derivations. More specifically, we develop a new operator algebra which involves some non- 
associative rules for the treatment of products of singular operators. Once this operator 
algebra has been carefully defined and developed, it provides a powerful tool for the formal 
manipulation of the equations and permits a careful and detailed comparison with the Bethe- 
Salpeter equations. It also alows us to derive several new results which would be difficult to 
derive using a purely diagrammatic approach. In applications the relativistic kernel for either 
the Bethe-Salpeter equation or the Gross equation is usually expanded in a perturbation 
series, and in this paper we obtain, for the first time, the form of the electromagnetic current 
operator for the Gross equation which is valid to all orders in this expansion. We also show 
explicitly that the theory conserves the charge of a bound state, and that gauge invariance 
is exactly preserved when the theory is truncated to any finite order, provided only that 
the strong kernel and the electromagnetic current operator are both truncated to the same 
finite order. 

This work is a continuation of recent work || in which the normalization condition 
for the three-body vertex function was derived, and also lays the foundation for extension 
of recent developments of the three-body Gross equations by Stadler and Gross JF]. The 
new algebra developed in this paper will be used to derive, in this forthcoming paper, 
the electromagnetic current operator for the three-body Gross equations |J, and we have 
developed the formalism here with an eye to this extension. Spectator currents have also 
been independently discussed by Kvinikhidze and Blankleider ||. Their discussion is more 
limited in scope than ours (here we develop an operator algebra, discuss the connection with 
the Bethe-Salpeter equation, and obtain results to all orders), but the results they do obtain 
agree with us (see the discussion in Sec. Ill below). 

A number of other works deriving the electromagnetic current for various relativistic 
equations have appeared recently. Coester and Riska have derived the current operator for 
the Blankenbecler-Sugar equation [JlO | and Devine and Wallace [if] and Phillips and Wallace 



12] have discussed the construction of a current operator for use with a relativistic version 



of the equal time equation. Extension of the new operator formalism presented here to these 
other equations is being studied. This effort may clarify a number of issues still unresolved 
in these treatments. 

This paper begins with a brief review of the Bethe-Salpeter equation and the correspond- 
ing current operator. In Sec. Ill we extend this discussion to the Gross equation, in both the 
unsymmetrized form for nonidentical particles and the symmetrized form appropriate for the 
description of identical particles. In Sec. IV we present the final form for the currents and 
show that the currents appropriate for identical and nonidentical particles are equivalent. 
We also show that the exact results in the two formalisms (BS and spectator) are identical 
if both are calculated to all orders. Then, in Sec. V we use the normalization conditions 
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FIG. 1. Diagrammatic representation of the integral equation for the four-point propagator. 

proved in a previous paper |J to show that the charge of the bound state is conserved 
by both theories. In Sec. VI we discuss the results when the perturbation expansions for 
the kernel and the current operator are truncated to a finite order, and show that gauge 
invariance is still satisfied. Finally, conclusions are presented in Sec. VII. 



II. TWO-BODY BETHE-SALPETER EQUATION 

In this section we review the Bethe-Salpeter formalism. Our results are not new, but 
the brief systematic development given here is needed both as an introduction to what will 
follow, and as a description of the formalism to which the spectator results will be compared. 
To prepare the way, we develop the subject using a conventional operator formalism. The 
need for non-associative operators will not appear until the next section. 

The operator form of the equation for the four-point propagator as represented in Fig. [I] 

is 

g = G BS -G BS vg (2.1) 

= G BS - GVG BS , (2.2) 

where the free two-body propagator G B s = —%G\Gi is defined in terms of the single-particle 
propagators Gi and V is the two-body Bethe-Salpeter kernel. 

The usual momentum-space forms of these expressions can be obtained by introducing 
the virtual momentum space states defined such that 

(x\ p) = e ip x ~ ip0t , (2.3) 

(p'\p) =(27r) 4 <*V-p) (2.4) 

and 

/(gjbMpl-l. (2.5) 

The operators are defined such that the momentum matrix elements for the one-body prop- 
agators are 

(ti\G i \p i ) = G i (p i )(27r)*5 4 (p f i -p i ), (2-6) 

the interaction kernel is 

(M\v\pm) = v(j/,n p)(2tt) 4 5 4 (p' - p) , (2.7) 
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FIG. 2. Diagrammatic representation the Bethe-Salpeter equation for the two-body scattering 
matrix. 



and the interacting two-body propagator is 

¥iP2\G\piP2) = G(p', P ; P)(2tt) V(P' - P) 



(2.8) 



where P = pi + P2 and P' = p[ + p' 2 are the total momenta in the initial and final states, 
and p = \{px — P2) and p' = \{p'i — p'2) are the corresponding relative momenta. 
The two-body propagator can also be written 



Q = Gbs — GbsMG 



where 



BS 



M = V - VGnsM = V- MGbsV 



(2.9) 



(2.10) 



is the two-body scattering matrix. The Bethe-Salpeter equation for the scattering matrix 
( |2.1(J| ) is represented by the Feynman diagrams of Fig. |2|. Equation ( |2.1| ) can be written as 



(G B l + v)g = i 

which implies that the solution for the inverse propagator is 

G~ X = G&+V. 

The equation for the Bethe-Salpeter bound-state vertex function is 

|D = -VG BS \Y) , 

which can be written 

= (1 + VG BS ) \T) = (G B l + V) G BS \T) . 



Using ( |2.12| ) this can be written 

Q- 1 l^) = o, 

where the Bethe-Salpeter bound-state wave function is defined as 

IV>> = Cbs |r> . 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



The scattering states are defined in terms of physical, on-shell states with the normal- 
ization 



(x\p) = e ip - x ~ iEpt 



(2.17) 
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FIG. 3. Feynman diagrams representing the five-point propagator. Inverse one-body propaga- 
tors are represented by the small, solid, square boxes inserted on the propagator lines. 



where E p = y/p 2 + m 2 . To include spin, we define the asymptotic single-particle plane wave 
momentum state as 



u {Pi s ) \p) f° r s P m 



\P) 



for spin = 



(2.18) 



The final state Bethe-Salpeter scattering wave function with incoming spherical wave bound- 
ary conditions is then 



■0 



(-) 



(Pi,si;p 2 ,s 2 \ (1 - MG BS ) 



(2-19) 



Using this 



g- 1 = ( Pl , Si; P2 , S2 \ (l - MG BS ) (G B l + V 

= ( Pl , Sl ; p 2 , s 2 | (G B l -M + V- MG bs V 



0, 



^2.20) 



where ( ^.10 ) and (pi, Si] p 2 , s 2 \ G B g = have been used in the last step. Similarly, the initial 
state scattering wave function with outgoing spherical wave boundary conditions 



[1 - GbsM) \pi,si,p 2 ,s 2 ) 



satisfies the wave equation 



Q' 1 V 



,(+) 



o. 



(2.21) 



(2.22) 



So the two-body Bethe-Salpeter wave functions for both bound and scattering states satisfy 
the equation 



Q- 1 



y-l 







(2.23) 



The two body current can be obtained from the five-point function describing the inter- 
action of a photon (with the photon leg amputated) with the interacting two-body system. 
This is represented by the diagrams of Fig. H, and corresponds to the operator equation 
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g* = -g (i jf g 2 - x + i j%g? + J£) ^ 



(2.24) 



where the inverse one-body propagators are introduced to allow for the factorization in terms 
of interacting four-point propagators. The inverse one-body propagators are represented by 
the square boxes inserted on the propagator lines in Fig. |3|. 
In order to demonstrate that the current 



J IA 



(2.25) 



is conserved, we must introduce the one- and two-body Ward-Takahashi identities in operator 
form 



q^Ji = \ei(q),Gi 1 



(2.26) 



and 



q^=[e 1 (q) + e 2 (q),V] (2.27) 

where ej(g) is the product of the charge (which might be an operator in isospin space) 
and a four-momentum shift operator defined such that 

(p'i\ei(q)\Pi) = e i (2 7 r)V( P : - Pi - q) . (2.28) 
Using the one- and two-body Ward-Takahashi identities give the following relation 



q^ = i fei^),^! 1 



e 2 (q),G 2 1 \ G^ 1 + [e 1 (q) + e 2 (q),V} 

i 



= e l {q)+e 2 {q),G B l + V\ = [ ei {q) + e 2 (q), Q 
This along with ( |2.23| ) implies that the two-body current is conserved 



J" 



0. 



(2.29) 



(2.30) 



For identical particles, the Bethe-Salpeter equation can be rewritten in an explicitly 
symmetrized form 



M — V — VGbsM = V - MG BS V 



(2.31) 



where V = A 2 V, M = A 2 M. and A 2 = 2 (1 + CP12) is the two-body symmetrization 
operator. (Note that Roman letters (e.g. M) are used for symmetrized quantities and script 
letters (e.g. M.) for unsymmetrized quantities, as in Ref. 0). The corresponding four-point 
propagator is 

G = A 2 G BS - G m VG = A 2 G BS - G BS MG BS (2.32) 

where G = A 2 g. The five-point function is also symmetrized in a similar fashion 

G M = -G (iJfGz 1 + iJ^G^ 1 + 7y G 

= -A 2 (G BS - GbsMGbs) (ufG? + UgG? + T ex ) (G BS - G BS MG BS ) (2.33) 

where J^ x = A 2 J£ X and satisfies the Ward-Takahashi identity 

q,T cx = [ ei (q) + e2(q),V}. (2.34) 

The proof of current conservation follows in exactly the same way as for the unsymmetrized 
case. 
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FIG. 4. The box diagram. 



III. THE TWO-BODY GROSS EQUATION 



In order to extend this discussion to the spectator or Gross equation, it is useful to 
examine the connection of the Gross equation to the Bethe-Salpeter equation. This is done 
most easily for the case of nonidentical particles. Identical particles will be discussed later. 



In order to introduce the singular operators needed for our discussion and to derive 
their non- associative operator algebra, we first review the procedure used to motivate the 
rearrangement of the multiple scattering series which leads to the Gross equation. This 
is illustrated by considering the second-order box diagram of Fig. f| which represents the 
interaction of two particles through the exchange of two light bosons. We assume the two 
particles to be of different masses, with the heavier mass associated with particle 1. The 
location of the 8 poles in the energy loop integration is shown in Fig. Here the positive 
and negative energy poles of interacting particles 1 and 2 are labeled 1 ± and 2 ± and the 
poles in the propagators of the exchanged bosons are unlabeled. For low energies the loop 
integral will be dominated by the the two poles 1 + and 2 + , which lie close to each other (and 
pinch above the scattering threshold). If the contour of integration is closed in the lower 
half-plane the result is dominated by the contribution from 1 + , the positive energy pole 
for particle 1. This suggests that it may be reasonable to separate the contour integration 
into two contributions, one containing only the contribution from the positive energy pole 
for particle 1 and one containing contributions from all of the remaining poles within the 
contour. 

This separation into two terms is best illustrated by considering the Dirac propagator 
for a single particle 



are the positive and negative energy projection operators. If we subtract the the Dirac 
conjugate 



A. Two-Body Equations for Distinguishable Particles 




(3.1) 



where 




(3.2) 
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FIG. 5. Location of the 8 propagator poles in the integrand of the box diagram in the complex 
Po plane (where po is the relative energy of the two internal particles). 
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FIG. 6. The singularities of the two contributions to the box diagram resulting from the decom- 
position of G\ into Q\ (left panel) and AGi (right panel). The role of the additional singularity 1 + 
in the upper half plane in the left panel is to pinch the contour. Mathematically this puts particle 
1 on-shell. 



A+(p) 



E P ~ P ° 
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from the first term on the right hand side of (|3.1|) and add it to the second term we obtain 



G 



d A p 



\p) 



m 



2zeA+(p) 



+ m 



E p (E p - p°Y + e 2 (E p -p° + le) (E p + p° - ie) 



(Pi 



(3.3) 



The first and second terms on the right hand side of this equation are represented by the left 
and right hand diagrams in Fig. |6], respectively. The first term contains a new pole which is 
the conjugate to 1 + , lies just above the real axis, and pinches the pole at 1 + when the limit 
e — > is taken. As we will see below, this automatically selects the positive energy pole for 
particle 1. The second term is a difference propagator corresponding to the second diagram 
in Fig. |^. It is the same as the original propagator but with pole 1 + moved above the real 
axis. If we now define 



Q 



IP) 



2eQ{p) 



(2vr) 4 2E p (E p -p°Y + e 2 



and 



AG 



d A p 



\p) 



Hp) 



(E p -p° + ie) (E p + pv 



(Pi 



(3.4) 



(3.5) 



where 
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N =S 2m for s P in = \ (3.6) 
1 1 for spin = , 




for spin 



i 



and 



i 



L(p)= U + m for spin- 
v ' II for spin = , 

we see that the propagator for particle i has been separated into two pieces 

G i = iQ i + AG i . (3.9) 

Furthermore, using contour integration it is easy to show that 

lim Q = l.m f *£ I ** \p) JL 2 j9^ (f) | 

-/<&£w«m«-?/<&£*«><m. (310) 

This shows that Q acts to place the propagating particle on mass shell and contains the 
projection operator Q = Q 2 on to positive energy spinor states, where appropriate. Be 
warned that Refs. || and did not make the distinction between Q and Q being made 
in this paper and that their Q is the same as our Q. However, because of the conventions 
( p. 17] ) to be introduced below (which were implicit in Refs. || and 0), this difference does 
not affect the conclusions previously reached in these papers and our results are consistent 
with these earlier references. 

While the introduction of the operator Q may seem straightforward, it is a singular 
operator and great care must be taken when using it. In particular, like the familiar delta 
function, its square is not defined. Later, we will be faced with the problem of how to treat 
quantities which naively appear to be products of singular operators, or a vanishing operator 
times a singular operator, and we will introduce a non-associative algebra for treating these 
products. Until then, the analysis is straightforward. 

Using ( |3.9| ), the Bethe-Salpeter equation (|2.10 ) for the t-matrix can now be formally 



separated into a pair of coupled equations. The first of these is 

M = U -UQ 1 G 2 M = U -UQ igi M, (3.11) 

or alternately 

M = U - MQigiU , (3.12) 

where U is called the quasipotential. The second equation relates the quasipotential to the 
BS kernel V. This equation is derived by requiring that the scattering matrix M. as given 
by ( |3.11| ) be identical to that of ( |2.1(J| ). The resulting equation for the quasipotential is then 
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FIG. 7. Feynman diagrams representing the Gross equation for the two-body scattering matrix. 
The cross on a propagator line designates that that propagator has been placed on its positive 
energy mass shell. 









e 






V 


+ 


V 




1 





FIG. 8. Feynman diagrams representing the quasipotential equation, 
propagator line represents the difference propagator. 



The open circle on a 



U = V - V(-iAG 1 G 2 )U = V - VAgiU = V - UA gi V . 
Note that we use the notation 

9i = G 2 
A 9l = -iAG 1 G 2 , 



(3.13) 



(3.14) 



where the propagator with particle 1 on shell is g\ = G 2 . [We find it convenient to label 
the two-body propagator by the on-shell particle and to distribute the singular factor of Qi 
which accompanies it to other parts of the equation (as discussed below). We have therefore 
introduced the lower case notation (i.e. g{) to distinguish the off-shell part of the two-body 
propagator from the one body propagator G 2 ] 

The pair of equations ( |3.11| ) and ( pM3| ), as represented in Figs. [7] and Fig. |, constitute a 
resummation of the multiple scattering series represented by ( |2.10| ) and are exactly equivalent 
to it by construction. The constrained propagator Q\gi in ( |3.11D limits the phase space 
available to particle 1 to the positive energy mass shell. Contributions from the remainder 
of phase space for particle 1 are included in the quasipotential (|3.13|) through the difference 
propagator Ag t . 

Using d3~T0|) in (|3~Tl| ) gives 



f d?p\ N 



(3.15) 



Note that the projector Qi has introduced a sum over all on-shell intermediate states for 
particle 1. In order to avoid the necessity of repeatedly writing the on-shell states and the 
associated sum, we will now introduce a notational convention. We will use the operator 
Qi to denote the presence of on-shell states acting on adjacent operators. If Qi appears 
between two other operators and therefore acts to both the left and right, on-shell states 
acting to both the left and right are assumed to be present. In addition the phase-space 
integral 
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is also assumed to be present. If Qi appears as the first or last in a string of operators and 
therefore acts to the right or left respectively, then only the corresponding on-shell states 
acting to the right or left are assumed. In this case no phase-space integral is assumed. That 
is, 

O'QiO^O'Y: I fl P 'Jl \ Pl ,s 1 )(p 1 ,s 1 \0 = 0'Q 1 
OQ 1 ^0\ Pl ,s 1 ) 

Q 1 0^( Pl ,s 1 \0 (3.17) 

where O' and O represent any nonsingular operators or string of operators. One consequence 
of this convention is the relation 

0'Q 1 Q 1 = O f Q 1 Q 1 = O f Q 1 0, (3.18) 



which follows from the observation that Q\ \pi, si) = \px, s\). Using the convention fl3.17| ) 
we can rewrite (|3.15[) as 



M = U -UQ igi M (3.19) 
= U-MQ l g 1 U. (3.20) 



We may also replace the Qi in Eqs. (|3.19|) and (|3.20|) by Qf; in this case the original Eq. (|3.11|) 



is recovered either by using the convention ( |3.17| ) on one of the factors of Qi and then using 
( 13.18| ), or, alternatively, by first replacing Q\ by Qi and then using the convention Q3.17|) . 



In Refs. and |7|] the Q used here was denoted by Q (and the conventions ( |3.17| ) and ( |3.18 
were implicit), so our results agree with those of these previous papers. 



Equation (|3.13|) represents a four- dimensional integral equation that is as difficult to 



solve as the original four-dimensional Bethe-Salpeter equation. However, as is shown in 
more detail below, this equation is usually approximated by iteration and truncation. Equa- 
tion Q3.19|) can be solved by noting that the constrained propagator Q\g\ requires that the 



scattering matrix on the right hand side of this equation has particle 1 constrained on shell. 
Replacing this using ( |3.20| ) gives 



M = U-UQ 1 g l Q 1 U + UQ 1 g 1 MQmU. (3.21) 

The fully-off-shell t matrix can therefore be obtained by quadrature from the t matrix with 
particle 1 constrained on shell in both initial and final states. This in turn can be obtained 
by placing particle 1 on-shell in the initial and final states in fl3.19| ) to give 

M n = U n - UngxMu 

= U n -M n g 1 U 11 (3.22) 

where Mn = QiMQi and U n = QxUQx. 
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In order to define the half-off-shell four-point propagator, we want to replace all of 
the propagators for particle 1 in ( |2.9|) with the on-shell projector Qigi. This can be done 
straightforwardly (i.e. avoiding the appearance of undefined factors of Qf) if the free particle 
inhomogeneous term is treated separately. We define 



'11 



Qi iG^ (G - G BS ) iG x l Q 1 + Q l9l = Q l9l - 9l M ll9l = Q l9l - g x U xx Q xx . (3.23) 



where the square brackets indicate that the propagators for particle 1 are first amputated 
using G, 1 and the result is then placed on shell. This equation for Gn can be written 

(9i~ X + Uu)Gn=Qx. (3.24) 

Since the projector Q x does not have an inverse, Gn does not have an inverse. However, 
the above expression indicates that Gn does have an inverse when acting on the subspace 
spanned by the physical particle states, ie. those projected out by the operator Q x . The 



solution of ( |3.24| ) on this subspace will therefore be written 

Gn =9i^ + U n , (3.25) 

where we bear in mind that Gn is defined only on the space spanned by the physical states 
of the first particle. 

The bound state vertex function for the Gross equation satisfies the equation 

|Ii) = -U n9l \T X ) . (3.26) 

This can be rewritten 

= (1 + U ngi ) \T X ) = (V 1 + Un) 9l |r x ) , (3.27) 

or 

GuWi) = 0, (3-28) 

where the Gross wave function is defined as 

|Vi) = 9i \Ti) . (3.29) 

The final state Gross scattering wave function with incoming spherical wave boundary 
conditions is defined as 

(tfj{' } \ = (pi,s 1 ;p 2 ,s 2 \(l-M n gi). (3.30) 

Using this 

'A~ } \ G xx = (Pi, sx;p 2 , s 2 | (1 - M n gi) (g^ 1 + U n ) 

= ( Pl , Sl ;p 2 , s 2 | (gi' 1 - Mu + U n - MugiUn) = (3.31) 



where ( |3.22| ) and (p x , si,p 2 , s 2 \ g\ 1 = have been used in the last step. Similarly, the initial 
state scattering wave function with outgoing spherical wave boundary conditions 
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FIG. 9. Box diagram with photon insertion. 



^! +) ) = (1 - 9iMn) \Pi, si; P2, s 2 ) 



(3.32) 



satisfies the wave equation 



Qn 1 lrf +) > = 



^3.33) 



So the two-body Gross wave functions for both bound and scattering states satisfy the 
equation 



Gn\i>i) = (i>i\Gn=0- 



(3.34) 



B. Two-Body Currents for Distinguishable Particles 

We now turn to the derivation of the two body current operator. This will be obtained 
from the five-point propagator as in our discussion of the BS equation. 

First consider the simple five-point box diagram shown in Fig. PI. The location of the 10 



poles in the energy loop integral is shown in Fig. |10| Since there are now two propagators 
for particle 1 in the loop, there are two positive energy poles for particle 1 labeled 1 + and 
1 /+ corresponding to the two propagators. If the contour is closed in the lower half-plane 
as shown in Fig. [II], the contour integral therefore contains two contributions corresponding 
to placing particle 1 on shell either before or after the photon absorption. The separation 
of propagators in the presence of the single nucleon current operator is then illustrated by 
the contour integrals displayed in Fig. |TT[ For spin 1/2 particles, the contour integral is 
decomposed into three terms 

00 dp Q f d 3 p /n ( m+ i>+\i \ TU , s ( m+ i>—\i 



J(^2K 0fA{p + l * q)J " ip+ > q) 




2 



d 3 p N ( m+ 



+ OfAdJfAdOi , (3.35) 

where Oi and Of are operators corresponding to the particle exchanges which occur before 
and after the interaction, p± = {E±,p) with E± = Jm 2 + (p ± hq) 2 , and the last term is 



the remainder of the dpo integration coming from all of the poles except 1 + and l' + . Note 
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FIG. 10. The 10 poles of the box diagram with photon insertion. 
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FIG. 11. Representation of the three terms resulting from the decomposition of the propagators 
of particle 1 in the presence of the one-body current insertion. In the limit e — > 0, the pinching 
poles in the top two figures insure that particle 1 is on-shell, either before or after the interaction. 
The bottom panel is the contribution from terms in which particle 1 is off-shell both before and 
after the interaction. 

that the singularities which appear in the first two terms after the integration cancel as 
q — > and that therefore the ie prescriptions have been dropped from the propagators. In 
algebraic form this decomposition can be written 

Of {dJ^Gx} Oi -> Of {Qi Jf Ad + AGiJf Qi + Ad Jf Ad} O t , (3.36) 

where the { } brackets indicate that only one loop integration is present even though there 
are two operators G\. 

Note that the expression (|3.36| ) does not contain the term QiJfQi which might be 
expected if the decomposition ( |3.9|) were blindly inserted into G\J^G\. In order to obtain 
such a term the contour integration would have to pick up the two poles at 1 + and 1 /+ 
simultaneously , which is clearly impossible. The only sense in which the contour integration 
might seem to pick up these two poles simultaneously is when they coalesce into a single 
double pole, which can occur for certain values of the external and internal loop momenta. 
However, even in these special cases the residue theorem 
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dz , = / </: 



C [Z - Z ) 2 JC 



Z - ZqY [z - Zq) 



2nif'(z ) 



shows that the only contribution comes from the single poles which result from the Laurent 
expansion of the integrand at the point Zq; there is no contribution from the double singu- 
larity itself. In our case, when the two poles do coalesce, the combination of the first two 
terms on the RHS (|3.36|) gives the correct result by producing a derivative term (similar to 
the f'{z ) term in the above example) arising from the cancellation of the singular parts of 
each term. 

Note that when the current couples to external lines, or when particle 1 is disconnected 
from the graph so that there is no loop integration involved, the term QiJfQi will be 
present. It vanishes only from internal loops. 

The relationship between various n-point functions as described in the Bethe-Salpeter 
formalism and the corresponding quantities for the Gross equation can always be obtained 
by a similar procedure. That is, starting with the Bethe-Salpeter quantity: 

1. Identify all loops contributing to the n-point function. 



2. Reduce all redundant products of one-body operators. For example in ( 2.24j) use 

3. In loops where the photon does not connect to particle 1, replace the one-body prop- 
agators for particle 1 with (|3] 



4. In loops where the photon does connect to particle 1, replace the quantity GiJfGi 
using ( |3.36| ). 

Careful application of this procedure will always result in a correct expression for the Gross 
n-point functions, and is straightforward when applied to the derivation of the Gross five- 
point propagator. However, in the application to three body systems it is necessary to treat 
the six- and seven-point functions, and the task of identifying all possible configurations of 
loops in these cases is quite tedious. In this case the task is greatly simplified if we develop a 
few identities which are equivalent to introducing a non-associative algebra for the operators 
which occur in the spectator theory. These identities also simplify the discussion of two- 
body systems, and will therefore be developed now. The discussion of the application of 
these ideas to three-body systems is postponed for forthcoming paper || . 

Since Q is very singular at the positive energy pole, considerable care must be taken in 
evaluating the product of this operator with other operators which may also be singular or 
vanishing at the pole position. To see this consider the product QiG~ l Q for scalar particles. 
Using O, 

lim QiG~ l Q = lim / \p') — ^ (p'\ iG~ l 

«-o e-o J 2rr 4 w 1 2E n , f RL/ - v' ) 2 + e 2 ^ ' 



(2tt) 4 2E p , (E p ,-p>°y + ei 

d 4 p , . _l 2e 

(2tt) 4 m 2E p (E p -p°) 2 + e* [Pl 
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lim 

e->0 



lim 



d 4 p 
d 4 p 



1 4e 2 



|P> 



(E p -p°Y + e 
4e 2 i (E p + p [ 



o 



it 



(2tt) 4 4£ p 2 (E p -jfl- ie ) (£ p _ p o + , e )2 
d 3 p iV 



(Pi 



IP) (Pi 



(2tt) 3 2E P 

and a similar result can be obtained for spin-| particles. This implies that 

iQiGrUQi -> 
A similar argument leads to the identities 

AGiG^AGi -> AG, . 



(3.37) 



(3.38) 



(3.39) 
(3.40) 



Note that these identities all refer to products where G i 1 is inserted between factors of 
Z\ = iQi or Zf = AGi, and can be sumarized by the compact statement 

Zf G ^zf^5 w Zf. 

However, repeating the derivation for operators O other than Zf gives new rules: 

a, 



AG i Gz 1 O i = O i Gr 1 AG i 
iQi G~ l Oi = Oi Gj 1 iQi - 



0. 



(3.41) 
(3.42) 



Hence AG, G^ 1 — > 1 and iQi G^ 1 — > for all operators except Z £ . These strange results 
can be understood if it is recognized that the operator algebra is not associative. When 
reducing products of operators the correct procedure is to first look for combinations of the 
form Zf iG^ 1 Zf and use rules ( |3.38| )-( |3.40| ) to reduce any which are present. After this 
is done, rules (|3.41|) and (|3.42|) can be used to further reduce the expression. Finally the 
conventions ( |3.17| ) can be used. These rules allow us to carry out formal operations on the 
operators which would be impossible or meaningless otherwise, and give us a truly algebraic 
way to obtain relations. As an example, using ( |3.38| ) permits us to show that the Qu given 
by the simple relation 



Si iG^giG^ Q^g 



-1 



11 



(3.43) 



is identical to that defined in Eq. fl3.23|) . The latter relation (|3.43|) has the advantage that 
it provides a more obvious and intuitive connection to the BS propagator. 



Finally, to implement the decomposition (|3.36| ) we introduce the rule 



Qi J? Qi 







(3.44) 



Using this, we can write 
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Gi J?Gi= {iQ x + Ad) Jf (zQi + Ad) 

-> zQi Jf Ad + Ad Jf Qa + Ad Jf Ad 

-> Qi ^ d + d ^ Qi + Ad Jf Ad (3.45) 



which reproduces ( |3.36|) . The rule ( 3.44 ) will always produce the correct result when used 



inside loops and when used to convert the combination zQiJfAd to zQiJfd in all 
connected diagrams. It agrees with the current derived diagrammatically by Riska and 
Gross and with the results obtained in Ref. ||. 

However, in a recent paper Kvinikhidze and Blankleider || have claimed that the last line 
in Eq. (|3.45 ) is in error, and they propose a "new" gauged propagator. Detailed examination 



of their result (see the Appendix A) shows that it agrees with the last line in Eq. ( |3.45| ) 



provided we treat the propagator d a $ a principle value, neglecting its imaginary part. But 
this is precisely the meaning of equations like ( J3.45 ). The role of the ie prescription in the 



propagator [e.g. as in 1/ (m 2 — p 2 — ie)} is to tell how to evaluate the contour integral over 
Po', once this integral has been evaluated and the result is no longer singular (which is the 
case for Eq. ( |3.45| ) where the singularities in each of the first two terms on the RHS cancel 



in the sum) we are instructed to set e to zero. In previously published work |J- || this 
was, in fact, done. Hence the results of Ref. |9j are identical to ours, and there is no error 
in Ref. @. 

We are now ready to use these new rules to reduce the BS five-point function with particle 
one on-shell. This is obtained from by first amputating the external factors of d and 
then placing particle one on-shell. This gives 

g& = Q x [iG^G* zdr 1 ] Qi = ~ (Q1G2 - QiG 2 MG BS ) (yf d 1 + 1J2G1 1 + J e M x ) 

x (QiG 2 - G BS MQiG 2 ) ■ (3.46) 

We now want to rearrange this expression so as to identify an effective current for use when 
particle one is on-shell. The basic procedure is to rewrite the five-point function so that it 
has a form similar to Q2.24j ), i.e. a four-point function with particle one on-shell, followed by 
a current, followed by a four-point function with particle one on-shell. This is accomplished 
by rewriting the above expression so as to include any contributions from the propagation of 
two off-shell particles within the effective current operator. To this end, consider the factor 

Qid - G BS MQiG 2 = Qid - (Qi - z'Ad) d-MQid 

= Qid - Q1G2MQ1G2 + iAG 1 G 2 (U - UQ 1 G 2 M) Q X G 2 

= (1 - A gi U) {Q igx - QigiMQm) (3.47) 

where in the first step the propagator for particle 1 is written in its separated form, and in 
the second step the scattering matrix is iterated using Q3.11J) . Similarly 

Qid - QidA^ds = (Q191 - QigiMQm) (1 - UA 9l ) . (3.48) 

This gives 

Glx = - (Qwi - Qi9iMQ l9l ) (1 - UA 9l ) (iJ^G, 1 + iJ^Gr 1 + J&) 
x (1 - A 9l U) (Q l9 i - Qi 9 iMQ l9l ) 
-> -{Qi9i ~ 9iM ngi ) Jfx (Q l9l - g x M n gx) = -G11J& On , (3.49) 
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where we used (|3.18|) in the last line and 

Jfi = Gi (1 - UA 9l ) (iJf G?a 1 + iJl'G'r 1 + J£ x ) (1 - A^C/) Qi (3.50) 
is the effective current for the Gross equation. This can be broken into two terms: 



j 1A,cS 



Q x (1 - UA 9l ) iJf Gj 1 + U^G?) (1 - A 9l U) Q 



■Ceff = Gi (1 " ^l) (1 " A^ET) fix 



(3.51) 



These forms will be used in our discussion of gauge invariance below. 

The effective current can be simplified. Using the rules ( |3.38| )-( |3.42p and ( |3.44| ) we obtain 



ii 



Qi 



J% - JfdU - UdJ? - iU(AG 1 J^AG 1 )G 2 U 



-lUAG^G^G^U + (1 - UA 9l ) J£ (1 - A 9l U) 



Qi 



(3.52) 



This form is convenient for calculations. 

As in the case of the four-point function (ft.llp , ( |3.49| ) is simply a resummation of the 
Bethe-Salpeter five point function (|2.33|) with particle 1 constrained on shell in the initial 
and final states. The two versions of the five-point function are equivalent by construction. 
This in turn guaranties that the matrix elements of the effective current between physical 
asymptotic states will also be equivalent. Any matrix element of this effective current is of 
the form 



J% - JfGxU - UG X J{ - iU (AGiJi AGi)G 2 U 
-tUAG 1 (G 2 J»G 2 )U + (1 - UA 9l ) J£ (1 - A 9l U) 



|V>i> • (3.53) 



We now show that the sum of the currents ([3.51 ) is gauge invariant. If we define 
U{G 2 1 + U^Gi 1 , Eq. (ET26D can be written 



% j ia = [ei(g) + e 2 (g),G B s] 



(3.54) 



Recalling ( 2.27 ) the divergences of the two parts of the effective current become 



^li.eff = Gi (1 - UA 9l ) [e x {q) + e 2 (q), G^] (1 - A 9l U) Q 1 
^eleff = Gi (1 - UA 9l ) [eM + e 2 (q), V] (1 - A 9l U) Q t . 

Adding these gives 

?m J u = Gi (1 - UA 9l ) [e^q) + e 2 (q), G' 1 ] (1 - A 9l U) Q x . 



(3.55) 



(3.56) 



Next we reduce the factor Q~ l (1 — A 9l U) Q\. The result we obtain depends on whether e x 
or e 2 multiplies from the left. If the factor is e\, use rules ( |3.41 ) and ( |3.42 ) and the equation 
for the quasipotential (|3.13j ) to obtain 
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e x g- 1 (1 - A 9l U) Q x = ei (iG^G^ [1 + i{AG x )G 2 U] + V - VA 9l u) Q x 

= e x {-U + U) Qi = 0. (3.57) 



A similar result holds for Q x (1 — UAgi) Q l , and we see that 



ii 



ei terms 







(3.58) 



independent of the fact that the initial and final states satisfy Eqs. ( |3.34| ). Hence Eq. ( |3.56p 
reduces to 



%J?i = Qi (1 - [e 2 (g), CT 1 ] (1 - A^) Qi • 



(3.59) 



To further reduce Eq. ( |3.59| ) we first use Eq. ( |3.13| ) to simplify terms involving the 
commutator [e 2 (g),y] 



^Jfi = Qi{[e 2 (q),G^} -UA 9l [e 2 (q),G^} - [e 2 (q),G^ s }A gi U + UAg l [e 2 (q),G B 1 s }Ag 1 U 

+ [e 2 (q),U] + U[e 2 (q),Ag 1 ]U^Q 1 

= Q 1 l [ tG^[e 2 (q),G 2 1 }+tUAG 1 [e 2 (q),G 2 }U+[e 2 (q),U} + U[e 2 (q),Ag 1 ]U^Q 1 

= Qi[e 2 (g),G 2 - 1 + f/]Q 1 = Qi[e 2 (g), Qn\Qi , (3.60) 

where the second equation was obtained using rule ( |3.39j ) to eliminate some of the terms 
linear in U and rule ( |3.40| ) to simplify the term involving Agi[e 2 (g), GggJAgi. The cancella- 
tion of the U 2 terms, leading to the third equation, then follows by substituting for Ag x and 
noting that AG X commutes with e 2 . Using the conventions ( |3.17| ), Eqs. ( |3.60| ) and ( p.34| ) 
imply that 



(3.61) 



so the current is conserved. 



C. Two-Body Equations for Identical Particles 

We will now extend the derivation of the two-body Gross equations to the case of iden- 
tical particles. Although simple arguments can be used to show that the result will have 
essentially the same form as those in the previous section with the substitution of appro- 
priately symmetrized quantities, we will proceed by considering a completely symmetrical 
approach to the construction of the four- and five-point functions. We will then show that 
necessary quantities can be reduced to a simpler non-symmetric form suitable for calcu- 
lation. In doing so we will illustrate the approach necessary for constricting the effective 
currents for the three-body Gross equation. 



Starting with ( 2.31 ) and making a symmetrical replacement of the one-body propagators 



in the intermediate state gives 
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M = V - V \ (Q igi + A gi + Q 2 g 2 + Ag 2 ) M 



(3.62) 



where g 2 = G\ is the propagator for particle 2 on shell, and Ag 2 = —iG\AG 2 . This can be 
rewritten as the pair of equations 

M = U - U \ (Q l9l + Q 2 g 2 ) M (3.63) 

and 

U = V-V\{/\ gi + /\g 2 )U. (3.64) 

There are now two channels that contribute the Gross equation, one where particle 1 is on 
shell and one where particle 2 is on shell. For the purposes of the following discussion it is 
convenient to pose the various equations in terms of a two-dimensional channel space. This 
can be done by introducing the vector 




(3.65) 



and the matrices 




(3.66) 



2 = ( t I ) (3.68) 



2 



and 



Q= Q nn)- (3-69) 



Q 



2 

We can now write 

Q l9l + A 9l + Q 2 g 2 + Ag 2 = D T (g Q + A^ ) D (3.70) 

The t-matrix equation is then 

M — U — UD T IgoQDM = U - UD T \g Q QD M , (3.71) 

where in the last step the limit e — > was taken, and the corresponding quasipotential 
equation is 

U = V -V '\D T 'Ag DU . (3.72) 
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Note that the factor 1/2 could be included in the definitions of go and Ago- We have chosen 
not to do this since the corresponding factors for the three-body case cannot be subsumed 
into the propagators. A closed form for the half-off-shell t-matrices is given by 

QDMD T Q = QDUD T Q - QDU \D T g QDMD T Q (3.73) 

Defining the t-matrix as a two-dimensional matrix in the channel space 

M = QDMD T Q (3.74) 

and the quasipotential in the channel space 

U = QDUD T Q , (3.75) 

the matrix form of the t-matrix equation is 

M = U -U \g M = U - M \g Q U (3.76) 

The nonlinear form of the t-matrix equation is 

M = U-M \g Q M - M \gJJ \g Q M . (3.77) 

Next the half-off-shell t matrix is parameterized in terms of a contribution from a bound 
state pole at P 2 = M 2 and a residual part 

ir) (ri 

where the bound state vertex functions are described by the vector of vertex functions with 
particle 1 or particle 2 on shell with 



|ri) 
|r 2 ) 



(3.79) 



Using the usual techniques, this gives the fully symmetrized two-body Gross equation for 
the bound state vertex function 



-U\g Q \T) (3.80) 



with normalization given by 



. , r ,(ldg 1 dUl \ 

1 = {Tl [2d^- 2 90 dP 2 2 90 ) lT) - (3 - 81) 

It is convenient to introduce the following definition for the interacting spectator prop- 
agator: 

9 = f </oQ - \9oM \g = \g Q - \g Ug = \g Q Q - gU \g . (3.82) 
This can be rewritten 
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(2<?o 1 + u)g = Q. (3.83) 
So the "inverse" of the propagator is 

g- 1 = 2g 1 + U. (3.84) 
The Gross equation for the bound state vertex function ( |3.80|) can therefore be rewritten 

= (l + E7±</ ) |r> = (20o 1 + U) I go \T) = g 1 |</>> (3.85) 



where the Gross bound state wave function is defined by 



|V>> = -so|r> = - 



l / l^i) 



(3.86) 



The final state Gross scattering wave function with incoming spherical wave boundary 
conditions is defined to be 



Using this 



(pi, si,p 2} s 2 \A 2 \D T (l - M\g Q 



1 



(-) 



4~ ] 



(3.87) 
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(-) 



(pi, si; p 2 , s 2 | A \D T (1 - M|sf ) (20Q- 1 + L/ 



(Pi, s a ; p 2 , s 2 \ A 2 \D T (20Q 1 - M + U - M\g u) = 



(3.88) 



where ( |3.76| ) and (pi, si; p 2 , s 2 | G" 1 = 0, for i = 1,2, have been used in the last step. Simi- 
lary, the initial state Gross scattering wave function with outgoing spherical wave boundary 
conditions 



satisfies the wave equation 



1 - ~0oM) D ~A 2 |pi, si; p 2 , s 2 ) = - 



(3.89) 



9 



V (+) ) = o 



(3.90) 



So the two-body Gross wave functions for both bound and scattering states satisfy the 
equation 



g-t\< l p) = (^\g-t=0. 



(3.91) 
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D. Two-Body Currents for Identical Particles 



Finally, we turn to the construction of the current for identical particles. Following 
the method previously developed, we obtain the current from the symmetrized five-point 
propagator for the Gross equation. This propagator is obtained from the symmetrized five- 
point propagator for the Bethe-Salpeter equation, Eq. ( |2.33| ), by replacing the two-body 
propagator, Gbs; associated with internal loops by the decomposition 



However, since the impulse term contains only one loop and the exchange term contains 
two loops, this substitution leads to a different result for these two cases. To illustrate this, 
consider the two-loop combination MGbs-Cx^bs^ which involves the exchange current. 
This combination gives 



where J^ x = DJ^D T . Note that the factors of 1/2 are the result of the fact that each of 
the two independent loops can be closed in two different ways. 

The comparable combination for the one-body current Gbs^ia^bs contains only one 
energy-momentum loop that can be closed in either of two ways so the symmetric separation 
of the propagators gives 



where Jj A = J IA 1. This argument shows that the factor J IA is transformed into 2Jj A . To 
complete the symmetrization we make the substitutions 



Gbs = \ (giQi + Agi +92Q2 + Ag 2 ) 



(3.92) 




(3.93) 



MG BS J? A G BS M = M\{( gi Q x + Ag x ) Jf A {g 1 Q 1 + A 9l ) 

+ (g 2 Q 2 + Ag 2 ) (g 2 Q 2 + Ag 2 )]M 
= M\D T (g Q + Ag ) Jf A (g Q + Ag ) DM 
= M\D T (g Q + Ag ) 2 Jf A ± (g Q + Ag ) DM 



(3.94) 



' 2 9oQ 



(3.95) 



for external two-body propagators, 



M 



DMD T 



(3.96) 



for the t matrix, and 




(3.97) 
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- -A 2 i flfoQ [l - £>M£> T | (g Q + Ag )] (2 Jf A + jy 

x 1-| (g Q + A 9o ) DMD T ] ±g Q. (3.98) 

As before, it is convenient to simplify the five-point function by incorporating any ap- 
pearence of off-shell two-body propagators within the effective current operator. To do this, 
consider the factor 

1 - \ (9oQ + A 9o ) DMD T = 1 - \g,QDMD T - \&g DMD T 

= 1 - \g Q QDMD T - ±Ag DUD T (l - \g Q QDMD T ) 
= (l - ±Ag DUD T ) (l - \g QDMD T ) . (3.99) 

Similarly, 

1 - DMD T \ (g Q + Ago) = (l - DMD T \g Q Q) (l - DUD T \Ag ) . 



(3.100) 



Equation (|3.98|) can then be rewritten 

G» - -A 2 |g Q (l - DMD T \g Q) (l - DUD T \Ag ) (2Jf A 
x (l - \Ag Q DUD T ) (l - \g,QDMD T ) \g Q Q. 

Using Eq. (|3.82 ) for the symmetric propagator, g, this becomes 

g^ = -A 2 gJ»g 

where the matrix current operator J 11 is given by 

J* = Q (l - DUD T \ Ag ) (2 Jf A + J^ x ) (l - § Ag DUD T ) Q 

— •* IA,cfT ^ J ex,eff ■ 



(3.101) 



(3.102) 



(3.103) 



We note for future reference that, using the rules ( |3.38| - |3T4"2"D and ( |3.44[ ), the contributions 
from the one-body current can be simplified, 



(3.104) 



(3.105) 



AgoJfAS 



AdJf 
AG 2 J 2 M 



Q 



(3.106) 



and 
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Ag Jf A Ag c 



Ad Jf AGiG 2 + G 2 J^G 2 AG 1 






AG 2 J^AG 2 G 1 + G 1 JiG 1 AG 2 



(3.107) 



We conclude this section with a discussion of the proof of gauge invariance for the 
symmetric current ( |3.103| ). The matrix form of one body Ward identity, Eq. ( |2.26| ), is 



<?m J ia = M?) + e 2 {q), G B l]l. 



(3.108) 



and using the two-body Ward identity, ( |2.27| ), together with the equation for the quasipo- 
tential (|3.72j) and rules (|3.38H3T4l|) , the four-divergences of the two parts of the effective 
current are 

q,J^ cS = Q2[e 1 (q) + e 2 (q),G B 1 s ]Q 

-Q[ei(q) + e 2 (q), Ag DUD T Q - Q DUD T Ag [e x {q) + e 2 (q), G B l} Q 
-QDUD T [e 1 {q)+e 2 {q) ) \Ag () }DUD T Q 1 (3.109) 

Q^es = Q[ei(q) + e 2 (q), DUD T ]Q+ QDUD T [ ei (q) + e 2 (q), \Ag Q }DUD T Q . (3.110) 

The terms quadratic in U cancel when the two equations are added, giving 

q^ = Q [ei(q) + e 2 (q), 2G B l + DUD T ] Q-Q [e x {q) + e 2 (q), G B l\ Ag DUD T Q 

- Q DUD T Ag [e\(q) + e 2 (q), G B g] Q . (3.111) 

This equation can be simplified using rules fl3.38| - j3l42"D 

q^ = Q[e(q), 2g Q - 1 + DUD T ]Q 



e{q),2g G - 1 Q + U = e(q),g 



where we have introduced the matrix charge operator, 

e(q) -■ 



e 2 (q) 
ei (g) 



(3.112) 



(3.113) 



Eq. ( p,112| ) is the symmetric generalization of Eq. ( |3.60|) , and leads, together with the wave 
Eq. (|3.91|), to the formal proof of gauge invariance. 



IV. CURRENT OPERATORS FOR THE GROSS EQUATION 

In the previous section we derived Eq. (|3.52|) of the current operator J{\, which is to 
be used for the treatment of nonidentical particles where particle 1 is on-shell, and the 
current operator J M , Eq. ( |3.103| ), for use with identical particles. In this section we will 
show that, using the symmetry of the states, ( |3.103j ) can be reduced to ( |3.52| ) [with the 
obvious requirement that the masses and charges are equal], so that the form (|3.52j ) can 
be used in both cases. We will then decompose ( |3.52| ) into individual terms and give a 
diagrammatic interpretation of the current. Finally, we compare our results with the Bethe- 
Salpeter equation. 
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FIG. 12. Feynman diagrams representing J^ t . The open circles on particle line 1 are the 
difference propagator AGi, the shaded rectangles are the quasipotential U, and the open rectangles 
with photon attached are J^ x . 



A. Equivalence of the currents 



First, we recall the simplifications of the symmetric one-body current terms given in 
Eqs. ( |3.104j - P".107| ). Using these results, the one body terms can be written in the following 
form 



J IA.eff 



Q 



2J£ + fx + jT 1 + j m jT + + f 



,4 



3% +3? +3" 



Q 



(4.1) 



where 



AP- - 

Ji - 

M 
Ji 

r-- 



iU 



AGi JfAGiGa + G 2 J^G 2 AG 1 + (1 <-> 2) 



U 



and j 2 is obtained from ji by substituting 2 for 1. Next, recall that (Vi 2 exchanges particles 
1 and 2 (where £ = ± depending on statistics of the particles), and use the identities 



|V>2> = CVl2 |Vl> 
A(? 2 = CPj2 &9l C^12 

Z7 = cp 12 17 = 17CPi 2 

Jl = cv 12 ji = r cx cv 12 = cv 12 r cx cv 12 



(V 12 U(V 



12 



to show that 
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FIG. 13. Feynman diagrams representing the matrix element of the effective current between 
bound states. 



Hence 



<-0i| jT|-0i> = (V'iIjTI^) = (HjZ^i) = (HtiM 
(A\ \A) = (Htf \i>i) = IV*) = 1^2) 



M r n 1^) + <Vii j? 2 IV2) + m j%x 1^1) + 



J22 1^2) 



UAG X G 2 J%G 2 U + (l - UA^) 7^ x (l - A^ZJ 



l^i) ■ 



(4.2) 



Note that this is identical to ( |3.53|) , provided that the symmetrized two-body interaction is 
substituted for U, the symmetrized interaction current is substituted for Jg x , and the masses 
and charges are set equal to each other. Hence we have show that Eq. (|3.53| ) may be used 
either for identical or nonidentical particles. 



B. Final expressions for the currents 



Now we will write explicit expressions for the matrix elements of the effective current 
between two bound states and between a bound initial state and a scattering final state. To 
facilitate this define an effective interaction current 



J? 



-iUAG l J^AG x G 2 U - iUAG l G 2 J^G 2 U + (1 - UAgA f ex (l - Ag x U 



(4.3) 



This current is illustrated diagrammatically in Fig. [12. 

The matrix element of the effective current between bound states can then be written 



Jr, 



(H J£ |V>i> - (H J?G 1 Ug 1 IT,) - (T 1 \g 1 UG 1 J? |^> + (^| |Vi> . (4.4) 
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FIG. 14. Feynman diagrams representing the matrix element of the effective current between 
a final scattering state and an initial bound state. 



In numerical calculations it is often convenient to introduce an off-shell vertex function 

|r> = -U gi IT,) , (4.5) 



where i = 1,2. This can be used to rewrite (|4.4j) as 

(VI J" IV) = <r| |Vi> + <-0i| Jfc?i |r> + (Vi| IVi) + <-0x] 4 l^i) • (4.6) 



The Feynman diagrams representing the elastic matrix element are shown in Fig. 13 . 

The matrix element of the effective current between a bound initial state and a scattering 
final state is 

(V H | J" IV) = ( P i,s 1 ;p 2 ,s 2 \A 2 (l - Mg 1 )Q 1 [jg - JfGjJ -UG X J>{ + j£ t ] |Vi) (4.7) 
Using the identities 



A 2 M = M 

A TV — JV- 
•^•2^int ~~ ^int 

(pi, s 1 ;p 2 ,s 2 \A 2 (1 - Mg x ) QiU = (pi, si,p 2 , s 2 \ M 
(Pi,si,p 2 ,s 2 \A 2 \ J% - JfGiC/J IVi) 

= ( Pl , Si; p 2 , s 2 | i[j^ 2 |r x > + jfd |r 2 ) + j£G 2 |r) + jfGx |r) 



(Pi,si;p 2 ,s 2 | 



JfGx |T 2 ) + j£G 2 |ri) =(pi,si;p2,s 2 | Jf IV2) + 



where to get the last relation we employed 



the current matrix element can be rewritten 
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V> ( = (pi,sf,p 2 ,s 2 



+Mg 1 Q 1 JZG 1 U + (1 - Mg 1 Q 1 )Jt nt |Vi> 



Pi, si;p 2 , s 2 Jf |V>2> + |V>i> - MG 2 J$ |-0i> - Jf |Vi> 



-MQ 1 J?G 1 G 2 \T) + (1 - MQiG 2 )J? nt 



(4.8) 



The Feynman diagrams representing the inelastic matrix element are shown in Fig. 14. 



C. Comparison to Bethe-Salpeter matrix elements 

Let us now compare the matrix elements derived above with those of the Bethe-Salpeter 
description. First, consider the elastic Bethe-Salpeter matrix element 



> BS = <r| g bs (u?g 2 1 + u%Gt + r cx ) g bs |r> 

= (r| - iG 1 J^G 1 G 2 - iG x G 2 J%G 2 - G l G 2 TjG x G 2 \T) 



(4.9) 



With the help of simple identities obtained from Eqs. (|3.9|) and (|3.44|) 



- %G X J{G X G 2 = Q 1 J[ U G 1 G 2 + G\J{Q\G 2 - iAG x J{AG x G 2 , 
—%G\G 2 J 2 G 2 = Q\G 2 J 2 G 2 — iAG x G 2 J 2 G 2 , 
-G 1 G 2 T CX G 1 G 2 = (Qi - iAG x )G 2 T ex G 2 {Q x - iAG x ) , 

relation Q4.5|) , and the conventions ( |3.17 ) and (|3.18|) , we can rewrite ( |4.9|) as 
M ^ I^)bs = (r| QiG 2 [J% - J?GJJ - UG X J? + J&] Q X G 2 \Y) 



(4.10) 



(4.11) 



with J(^ t defined by ( |4.3| ). We have demonstrated again that our spectator matrix element 
Q4.4D exactly equals Bethe-Salpeter one ( |4.9| ). Of course, this should be so by construction. 
Still, the derivation of this section gives a useful shortcut to the correct spectator matrix 
element. It also illustrates how the current J 2 — J X G X U — UG x Ji in 



and ( 4.11 ) follows 

from the Bethe-Salpeter matrix element if one puts the first particle on-shell, i.e., if one 
keeps only the first term in Eq. ( [3.72 ) for the quasipotential (for consistency we should 
replace U — > V, which also holds when all terms with AGi are omitted). The last part of 
the effective current J£, t then gathers all higher order effects. A very similar consideration 
can be applied to the break-up matrix element (|4.7|) . One finds that the parts of the matrix 
element with loops can again be related by identities ( |4.10| ), while the loopless parts, i.e., 
an IA contributions without a final-state interaction, are identical for both approaches. 
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V. CHARGE CONSERVATION 



In this section we show that the total charge of a bound state is equal to the sum of the 
charges of its constituents, e\ + e 2 , and discuss how this result emerges automatically in the 
Bethe-Salpeter and spectator formalisms. In this discussion we will assume for defmiteness 
that the two particles are nonidentical, but our results will hold for identical particles also 
since the current operator in the latter case is obtained by symmetrization of the former 
one. 

First recall that taking the q — > limit of the one-body Ward-Takahashi (WT) identity, 
Eq. ( |2.26|) , implies that the one-body currents satisfy 



Jt(0) = ~e 



dpiu 



(5.1) 



This relation will be used in both formalisms. 

Next consider the two-body WT identity, Eq. (|2.27 ), in the context of the BS formalism. 
It is well known [I3| that the contribution to the charge operator which comes from the 
exchange current can be uniquely determined by taking the q — > limit of (|2.27|) . The 
derivation of this result was discussed in great detail by Bentz [14], and we only review 



it briefly here. Since the overall four-momentum is conserved, the kernel is a function of 
only three independent four-momenta, which can be chosen to be either P, p, and p', or 
P, pi, and p' x . Depending on how we choose the independent momenta, the q — > limit of 
Eq. ( PT27D gives 



^x(O) 



ei + e 2 



-(ei + e 2 ] 



dV(p',p,P) (ei-e 2 ) 

dP, 2 
dVjp'^P) 
dP u 



dV{p',p,P) dV{p',p,P) 



dV^px.P) 



dp[ 



+ 



dPf, 
dVjp'^puP) 
dpi,„ 



(5.2) 
(5.3) 



where in (|5.2|) the partial derivative with respect P M implies that the independent vectors p^ 
and p'^ are held constant, while in (|Q| ) the partial derivative with respect P M implies that 
the independent vectors p\^ and p[ are held constant. Similarly, 



— 1- 



.dGT l ( Pl ) 



dp 



dG- 1 



13S 



(p,p) 



IdGvl 



(P,P) 



dP^ 
dGgj^P) 
dP„ 



2 dpp 
9Gzk(PuP) 



dm 



and 



iG^ipt 



dGn 1 



{P2 



dp 



dG£(p,P) + ldG£(p,P) 



dP» 
OP, 



dp. 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



The correct forms of these equations depend on our choice of independent vectors. In the 
BS case either of the forms can be used since there are no additional constraints on the 
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vectors, but in the spectator case with particle 1 on-shell we must use Q5.3Q , Q5.5Q , and 
because p 2 will explicitly depend on P in cases when pi is constrained. 

We use ( p.2| ), (|5T5|), and ( p.7| ) to evaluate the bound state matrix element of the charge 
operator in the BS formalism 

(V>| no) m = - m ex i^p^G- 2 \ P2 ) + e 2 iG^\p x )^p^ - J£(0) \^) 

(dGv\(p,P) , dV(p',p,PY 



[ei + e 2 ] 



V dP^ dP^ 



| e 1 -e 2 / dG B 1 s (p,P) + dV(p',p,P) + dV(p',p,P) 
2 V dp^ dp'^ dp^ 

{e l + e 2 )2P» 1 (5.8) 



where the normalization condition for the BS vertex function P,|13|l was used in the last 
step to simplify the (ei + e 2 ) terms, and the cancellation of the (ei — e 2 ) terms follows from 
integrating dG^/dp^ by parts and using the bound state BS equation ( |2.23| ). The final 
form of (|5.8| ) shows that the charge is conserved. 

Now we turn to the spectator formalism with the effective current given in Eq. fl3.52| ). 
We begin by pointing out that, unlike in the Bethe-Salpeter case, one cannot obtain Jn(0) 
from the corresponding Ward-Takahashi identity ( |3.60| ). For nonidentical particles the clear 
indication of this fact is that the charge of the first particle (which can be completely 
arbitrary) is absent from the WT relation (|3.60|) and any current determined from this 



relation would therefore depend on e 2 only, which is certainly not correct. The reason for 
this was alluded to in Sec. Ill: the condition restricting the first particle to its mass shell 
leads to an effective current in which the terms proportional to the charge of particle 1 are 
purely transverse. There can be also transverse currents in the BS case, but they are of the 
form a^q", with antisymmetric and nonsingular for q —* 0, and hence they vanish in 
this limit, and all parts of the current contributing to the charge can be recovered from the 
WT identity (see Ref. [0). In the spectator formalism those parts of the current which are 
transverse by virtue of the on-shell condition do not vanish in the q — > limit. Therefore, 
the effective current in q — > limit cannot be fully recovered from the WT identity and has 
to be obtained by taking the limit explicitly. 

The effective spectator current for zero photon momentum follows from Eqs. (|3.52|) and 

( dG^ 1 

^i(O) = 2i{- e 2^; - i J i G i u + UG 1 Jn q ^ + (1 - UA 9l ) J£(0) (1 - A 9l U) 



+ 



dG^ 1 dG^ 1 
iei uAGi—1— AG 1 G 2 U + ie 2 UAG l G 2 —^- G 2 U 



Qi- (5.9) 



The first term in the last line of this equation can be reduced if we use rules ( |3.40| ) and 
( |3.41| ) and integrate by parts twice (noting that p\ is unconstrained in this loop and that P 
is to be held constant) 

ie.UAG^ A Gl G 2 U = -i ^f^ AG^V - ie^AG ^f^^ 
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-2ie 1 U^-G 2 U-ie 1 UAG 1 ^-U 
e 1 U^U + ie l UAG 1 ^U, 



(5.10) 



where here and below p[ and p x are the four-momenta of particle 1 after and before the 
interaction, respectively, and p'[ denotes the momenta of the loop integration implied by the 
product U . . .U . Using 



dG 2 \p 2 ) dG 2 (P- Pl ) 
<-J2 ; ^2 



dp2,n ' dp 1;ll 
the second term in the last line of Eq. (|5.9|) becomes 



te 2 UAG 1 G 2 ^f- G 2 U = ie 2 UAG 1 ^U . 

dp 2 ^ "V\,n 

The term with the exchange current J^(0) is simplified with the help of 

dU dV 



dp'i^ 
dU 
dp hll 



dp[ 



[1-A gi U) 



[1-UAgS 



dV 



dp 



(5.11) 

(5.12) 

(5.13) 
(5.14) 



These relations are obtained by differentiating the corresponding off-shell quasipotential 
equations and using the fact that the structure of the integral equations insures that the 
only dependence on the final momentum p\ in ( |5.13|) , or on the initial momentum p\ in 
( f5.14j ), is found in the kernel V. A similar argument gives 



OU _dV ri _ AgiU) _ vd (A 9l U) 



dP, dP^ 



dP„ 



(5.15) 



and hence 



/ TTA . dV „ . rn dU TTA dU TT d(A gi U) 
1 - UA gi ) — (1 - A 9l U) = — - UAg x — + U- K 



dP„ 



dP^ ^dP„ 

dU + u d -^u 



dPa 



BP, 



BP, 



(5.16) 



Hence using (p|), fl5l3p , flOj) , and (1516]) gives 



1 - UAg x ) J£(0) (1 - Ag x U) 

dU 

-(ei + e 2 
+ei 



<9R 



ei 



UA 9l A 9l U 



dU(p' l ,p 1 ,P) + dU{p' 1 ,p 1 ,P) 



dp K 



(ei + e 2 



t7 — (7 — ei 



3P, 



+ T! h U^r-H-U 



dp'! 



(5.17) 
(5.18) 
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Since p" is off-shell in the integration loop, we could integrate by parts to simplify the last 
line of ( gT7l) . 

Making these substitutions and combining terms permits us to simplify Eq. (|5.9|) 



^i(O) = Qi{-e 2 ^- - [J^G 1 U + UG 1 J^ 



- (ei + e 2 



dU 
OP, 



( du( P [, Pl ,p) du(p' llPll py 



' 1 I dp' 1>fl + dp 1>fi 
dG 2 dA 9l 



+ (ei + e 2 ) U 



iAG v 



U Q 



dp'i, dP, J 

Recall that partial derivative with respect to P holds p\ constant, and hence 
dA 9l ( Pl ,P) _ ^ AGM dG 2 {P-K) = iAGi{pi] 9G 2 (P- Pl ) 



(5.19) 



BP, 



dP„ 



9 Pi„ 



(5.20) 



so the last line proportional to U 2 in (5.19) cancels. The remaining terms will be now shown 
to be proportional to the normalization condition. Let us point out that exactly such terms 
(with U — > V) would appear if only the leading order quasipotential with corresponding 
currents are considered. 

To simplify (|5.19|) one has to reduce the term [J±GiU + UGxJ\] . It must be treated 
carefully as each term is singular as q — > 0, but, as discussed in Ref ||, the singularities 
cancel in the sum. Following the argument developed in Ref. || , using the notation p\ and 
p' x to indicate those cases where the four-momenta of particle 1 are restricted to their mass 
shell, and exploiting the bound state equation (|3.28| ) gives 

- (^1 G 2 [J*£G X U + UG x J?\ q _ Q G 2 |r x ) 

= <r x | G 2 [UuGztfdU + UG x J?G 2 U n \ G 2 \Tx) q _ 

* -- (TilGJUufafa P + q) Wipl) ~ G^ip'l - q)] G x ip'{ - q) 

xG 2 (P - p'i + q)U{p'{ - q,px, P) + U(p[,p'l + q,P + q)G 2 (P - p'[) 

xGxip'i + q) [G^ip'i + q) - G:\p1)} Uxxip'lpx, P)]g 2 \T x ) q ^ 

= — <r x | GJUnip'xJl P + q)G 2 (P - p'i + q)U{p'[ - q,p 1 , P) 

-U{p[,p > ; + q,P + q)G 2 {P-p^U 1 x(plp 1 ,P)}G 2 \T 1 ) q ^ 



-G 2 Un — U\\G 2 )G 2 \Ti) 
dpi^' 



ei (T, 



dG 2 (P-px) 
OP, 



-G 2 (P - p[) ( d -m^H + dU ^ P ) \ g 2{P _ Pl) lTl) , (5.21) 
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where, in going from the second to the third step we used G^ l Qi = for the on-shell 
momentum p", and in the last step we used the bound state equation to remove the factors 
of G 2 U whenever possible. When simplifying (|5.21 ) it is important to choose the dummy 



integration momentum p'[ so that the on-shell condition does not depend on the photon 
momentum q. Finally, substituting this result into ( |5.19| ) gives 



•tfi(O) = (ex + e 2 ) Q l \g?^G? - |^| Q l , (5.22) 



and the elastic matrix element of the effective spectator current at q = becomes 

Jfi(o) |Vx) = <r x | G 2 jf 1 (o)G 2 |rx) 

'dG 2 (P- Pl ) 



(ei + e 2 ) (r\ 



-G 2 (P - £) G 2 (P - ft) ) |rx) . (5.23) 

However, the normalization condition for the spectator vertex function is just 

<r, I ( 9G2 <f p ; f " ) - C 2 (P - 9t/( ^- P) G 2 (P - ft) ) |r,> = 2P» . (5.24) 

This was discussed in great detail in 0, where it was derived from the nonlinear form of 
the spectator equation without reference to the e.m. current (in that reference the spectator 
kernel was denoted by V, but the derivation did not specify the kernel in any way and holds 
equally well for the kernel U). Obviously the relations ( |5.23| ) and ( |5.24 ) are consistent with 



<r x | G 2 J^{0)G 2 |r x > = (e a + e 2 ) 2P» , (5.25) 

which is the statement that the charge of the bound state is ex + e 2 , completing our proof. 

Our derivation is valid for any interaction V and the corresponding quasipotential U (e.g., 
also for phenomenological ones, such as a separable interaction). It is only necessary to have 
an interaction current at the Bethe-Salpeter level consistent with the one body current, so 
that the total BS current is conserved. Furthermore, since we have not specified the spins 
of the constituents or the bound state in our derivation, it should apply for arbitrary spins. 



VI. TRUNCATION 

To this point, no approximations have been made in constructing either the n-point 
functions or the effective current operators. In particular, the equivalence between the 
Gross equation and its quasipotential and the Bethe-Salpeter equation is exact only if the 
Bethe-Salpeter kernel V and the spectator kernel U are related by Eq. ( |3.13| ). This means 
that if one of the kernels is truncated to some finite order, the other must involve terms 
of all orders. In practice, both kernels are generally truncated to some finite order and the 
two formalisms do not give identical results. The usual approximation is to keep only the 
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one-boson-exchange-contribution, either for or [A 1 ). The problem is then to verify that 
the various relations leading to conserved current matrix elements are maintained in the 
presence of the truncation. 

First assume that we have some Bethe-Salpeter kernel V — > XV and the associated 
current J M — > Jf\ + A where the interaction current and the interaction satisfy the Ward- 
Takahashi identity ( |2.27| ). [In this section we will again limit the discussion to nonidentical 
particles.] Here the parameter A has been introduced to assist in the counting of occurrences 
of the interaction V and the associated exchange current J^ x and will eventually be set to 
unity in all calculations. From ( |3.13| ) it is clear that the quasipotential can be written as a 
series in A as 

oo 

[/ = ]T X N UW (6.1) 

N=l 

Substituting this into ( |3.13j ) and collecting the coefficients of the various powers of A, we 
can identify the quasipotential of the iV-th rank U^ N ' as 

U {1) =V } (6.2) 
UW = -VA gi U( N -V , N > 1 . (6.3) 

Similarly, using (|3.50|) implies that the effective current can also be expanded 

oo 

J& = E A^jg " (6.4) 

N=l 

where the iV-th rank contributions to the effective current are .7^ ^ = J^veff + ^exls- For 
JiA eft t nese contributions are 

^& = Qi^Qi, (6-5) 
4Zs = ~Qi (U {1) A 91 J? A + J&AgtUV) fix , (6.6) 



/ iV— 1 \ 

= ~Qi \U {N) ^9iJh + JU^ N) - E U^Ag^Ag^ Qx , 

V M=l / 

if N > 1 , (6.7) 

and 

4% = 0, (6.8) 

r(l)M 



= QiJ&i, (6-9) 
4x!es = -Qi (u {1) a 9i j^ + J? x A9iU {1) ) Qi, (6.10) 

& = -fii [U^A 9l J^ + J^ x A9iU^ - E U^^Ag^Ag^A Qx , 

if N>2. (6.11) 
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At the lowest rank N = the particles do not interact and only disconnected diagrams 



[which are not fully described by the current ( 3.50 )1 occur. To get a nontrivial description of 
interacting particles and their effective currents one has to include at least the rank N = 1 
terms. 

It is easy to show that a theory truncated at rank N is gauge invariant (and also covariant 
of course) provided all terms up to and including rank N are included. To do this we use 
Eqs. CT-QB^D and rules (ggg)-(ggg) in Appendix B to show that 



(0)M 

11 



+ 4?) 



Qi [e2(q),G^ 1 + uS ) ]Q 1 
Qi e 2 (g),f/ 1 ( f ) l Qi 



(6.12) 
(6.13) 



where ( |6.12| ) holds for the sum of A = and A = 1 terms, and ( |6.13 ) for any finite 
A > 2. Hence all terms linear in e\ cancel exactly in the truncated WT identities, (|6.12 ) 
and ( |6.13| ), just as they do in the untruncated identity, Eq. ( |3.6U| ), and the results ( |6.12D 
and ( |6.13| ) are completely consistent with ( |3.60| ). We have shown that an effective current 
which is the sum of terms up to any rank A max > 1 is gauge invariant provided only that the 
quasipotential and the current include all contributions up to rank A max . Furthermore, the 
derivation required only that the BS kernel and BS current satisfy (|2.27| ); they are otherwise 
unspecified. 

Now consider a Bethe-Salpeter potential consisting of two independent contributions 



V = \ 1 V 1 + X 2 V 2 



(6.14) 



with corresponding exchange currents AiJ^ ex + \ 2 J 2<eyi where the two components of this 
current satisfy (|2.27|) with the corresponding components of the potential. Examination of 
( j3.13j ) indicates that the quasipotential can be expanded in the form 



jVi,iV2=0 



where Eq. Q3.13| ) gives 



= v 2 , 

t/W.o) = -VtAgtU^- 1 ^, 
jj(o,n 2 ) _ -V 2 Ag 1 U^°' N2 ~ 1) , 



N t > 1, 
N 2 > 1, 

-V 2 A 9l U^ N ^ 



Nx,N 2 > 1 

Using ( 3.51| ), the corresponding contributions to the effective current are 



t(0,0)m _ n T n n 
■Jih ff — feil^iA ^1 



J, 



IA.eff 

(N u N 2 )ii 
IA.eff 



AT 2 ) 



(6.15) 



(6.16) 
(6.17) 
(6.18) 
(6.19) 
(6.20) 
(6.21) 



(6.22) 
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Ni N 2 . 

£ £ U( N i- M i> N *- M ^A gi J? A A gi U^> M >nQi , N l0 i N 2 >1, (6.23) 

Mi=0Af 2 =0 ' 

-&t = , (6-24) 
•&=QiJ?<*Qi, (6-25) 
j£$? = QiJZ-Qu (6-26) 

+ J^A^t/^- 1 -^) + J» ex A gi U^ N *-V 

Ni-l N 2 

-EE U^ Nl - Ml - 1 ' N9 - M ^Ag l J^Ag 1 U^ M ^ 

Mi=0 M 2 =0 

Ni Na-1 } 

-E E U {Nl - Ml ^- M '- 1) Ag 1 J^ cx Ag 1 U^ Ml ' M ^ \Q U N 1 ,N 2 >1. (6.27) 

Mi=0Af 2 =0 J 

The divergence of the effective current is then 

q,j[i 0) " = Qi [e2(q),G?] Q u 
qAi M = Qi[zMMi 1 ' N2) ]Qi N i0 tN 2 >1. (6.28) 

This implies that if all terms up to A^i max and iV2 max are retained in the quasipotential and 
the effective current that the Ward-Takahashi identity will be satisfied. Note that iV lmax and 
-^2max do not have to be equal. That is, contributions from the two parts of the interaction 
can be truncated at different orders without disturbing the Ward-Takahashi identity. 

The implication of these two results is that it is possible to truncate the quasipotential 
and interaction current in a consistent fashion without disturbing the Ward-Takahashi iden- 
tities and that the truncation can happen at arbitrary orders. Indeed, from this it is clear 
that the requirement of current conservation places little constraint on the truncation of the 
equation. Some other physical consideration must then determine the method of truncation 
of these quantities. 

An often used approximation to the Bethe-Salpeter equation is to collect contributions to 
the kernel containing the same number of boson exchanges. This is a natural procedure in the 
case of a perturbative approximation for a weak coupled field theory. This approximation is 
also used in relativistic models of the nucleon-nucleon system where the justification is that 
irreducible contributions with increasing numbers of exchanged bosons have a shorter range 
and tend to have a small effect on the wave functions and low energy scattering amplitudes. 

Consider an interaction following from multiple exchanges of the single type of boson 

oo 

V = ]T V {n) , (6.29) 

n=l 

where the superscript n denotes the number of exchanged bosons and is an irreducible 
contribution to the Bethe-Salpeter kernel. Again, from ( |2.27|) it follows that the Bethe- 



Salpeter exchange currents can be decomposed in a similar way 



37 



7M _ T(")M 
"ex / j °ex ' 



(6.30) 



n=l 



and the Ward-Takahashi identity is satisfied separately for each n. Actually, in passing to 
our quasipotential framework we can formally consider each set of Bethe-Salpeter-irreducible 
contributions of fixed n to be independent contributions in the sense considered in the second 
case discussed above. The quasipotential and effective current for each contribution could 
then be truncated independently of the others. 

However, it has been shown that the convergence of the Gross equation is improved, in 
some cases, by a delicate cancellation of crossed-box diagrams and subtracted box diagrams 
of the same order in n arising from the iteration of the quasipotential equation. Therefore, 
the physical consideration of convergence may require that contributions to the quasipo- 
tential with a fixed number of boson exchanges also be collected together. That is, the 
quasipotential can also be expanded 



U = ]T U {n) 



(6.31) 



n=l 



where n is the number of exchanged bosons contributing to U^ n \ Substituting this into 
(EM gives 



[/(i) = yd). 



n-1 



Jj{n) = y(n) _ £ y(n-a) A 9l U {a \ 71 > 1. 



o=l 



Using (|3.51| ), the corresponding contributions to the effective current are 



IA.eff 



(6.32) 

(6.33) 



(6.34) 
(6.35) 



IA,efT 



Qi {UVA 91 J& + J&AgxUW) Qi , 
-Q l ^Ag l Jf A + J? A Ag l U^- n ^ n>l, (6.36) 



J, 



(0) M 



cx,cfT 



J. 



ex,eff 



J. 



(2)/i 



ex,efT 



J, 



(n)n 



ex,efF 



( n—1 n— 1 



'A ■ 



(6.37) 
(6.38) 
(6.39) 



a=l 
n—2 n—a—1 



a=l 



n > 2 . 

a=l 6=1 J 

The divergence of this effective current is (see Appendix B) 

Qi fe 2 (g),G 2 - 1 l Qi 



Hli J ll 



e 2 (q)Mi 



in) 



n > 1 . 



(6.40) 



(6.41) 
(6.42) 
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This implies that the Ward-Takahashi identity is satisfied if the quasipotential and effective 
current include all contributions from boson exchanges up to some n max . This can be easily 
generalized to include additional kinds of bosons. From the second case presented above it 
is also clear that the equations can be truncated at different numbers of boson exchanges for 
each type of boson. For example, a meson exchange model of the nucleon-nucleon interaction 
could contain contributions from up to two pion exchanges, but heavier meson contributions 
could be truncated at the one-boson-exchange level. 

VII. CONCLUSIONS 

This paper develops a detailed algebraic treatment of the spectator or Gross description 
of strongly interacting two-particle systems in the presence of an external electromagnetic 
field (treated to first order). Our factorization of the five-point function follows naturally 
from the original definition of the spectator equations. 

We start from the Bethe-Salpeter formulation, i.e., we assume that the underlying dy- 
namics is known in principle and that it generates a series of Feynman diagrams which 
specifies both the interactions of two-nucleon system (Bethe-Salpeter equation) and the in- 
teraction of the two-nucleon system with an external electromagnetic field (Bethe-Salpeter 
exchange currents). The Bethe-Salpeter currents satisfy a Ward-Takahashi identity involving 
the Bethe-Salpeter four-point propagator. 

The spectator description is shown to result from rearranging these sets of diagrams, 
expressing the dynamics effectively in terms of a modified free two-nucleon propagator: in 
intermediate states one of the nucleons is restricted to its positive energy mass shell. The 
parts of the original diagrams in which this constraint does not hold are summed into a new 
effective interaction kernel (quasipotential) and an effective current (interaction current). 
The effective current satisfies a Ward-Takahashi identity with the corresponding four-point 
spectator propagator, so that the current is conserved. When all terms are included, the 
wave functions and current matrix elements are identical to those of the Bethe-Salpeter 
formalism. 

In applications, the whole infinite set of diagrams is not generally included, and we show 
that the series can be truncated to any finite order and still preserve gauge invariance. 
Most applications of the Gross formalism have been made using the lowest (second-order) 
one-meson-exchange approximation. Formally, this paper defines a consistent formulation 
for any finite order, and also shows that it is possible, for example, to include consistently 
the forth-order two-meson exchange contributions for some of the more important mesons 
(perhaps only the pion) while at the same time limiting the treatment of heavier mesons to 
the lowest, second-order. 

Although we have confined the arguments of this paper to the construction of electro- 
magnetic current matrix elements, the method is general and can be used, for example, 
to treat weak and axial vector currents. The extension of this formalism to three-particle 
systems will be presented in a future paper ||. 
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APPENDIX A: THE ONE-BODY CURRENT FOR PARTICLE 1 

In this Appendix we briefly discuss the comments of Kvinikhidze and Blankleider || in 
more detail. 

To illustrate the issue, consider the following contour integral 

^ 1 f f(z) dz , , , 

1= — /- J -\r? r, A.l 

2ixi Jc (z — z\ — iei)(z — z 2 — i&i) 

where the contour C encloses the two poles at z\ and z 2 , f(z) is analytic inside of the 
contour, and the limit — > is implied. Evaluation of the integral is straightforward, and 
gives 

1 = ^T-r{M + - + ie 2 )} 

Z\ - z 2 + ibe L J 

-> /'(zi + ie x ) -> /'(^) as z a -> 22 , (A.2) 

where the contour C encloses the two positive energy poles only, be = e\ — e 2 . Note that 
zero in the denominator at z\ — z 2 + ibe = is canceled exactly by a zero in the numerator, 
so the final result has no singularity. 

In the derivation of the one-body current for particle one, leading to Eq. (|3.45| ), we are 
confronted with a similar integral. In that case, in the Breit frame, the integral comparable 
to ( |PD is 



jp = j dko Jx{k + \q,k-\q) 



c 2ni (E + + k - iex)(E_ + k - ie 2 )(E + — k — ie 1 )(E_ — k — ie 2 ) 

1 f J?(K + jg, K - \q) Jf(A;_ + \q, k_ - \q) 

( E _ - %e 2 f - (E + - ie 2 ) 2 \ 2(E + - ie x ) 2(£_ - ie 2 ) 



(A.3) 



where the contour C encloses the two positive energy poles only, E± = ym 2 + (k ± |q) 2 , 
k + = {E + — iei,k), and A;_ = (E_ — ie 2) \t). Once again, the zero in the denominator at 
E_ — E + + % be = is canceled exactly by a zero in the numerator, so the final result has no 
singularity. However, the first two terms in the last line of Eq. Q3.45| ) (identical to the last 
two terms of Eq. (2.33) in Ref. [BJ), in the notation of Eq. ( |A.3| ), become 

O J^C +C /"O ~ f Ji( k + + h k +~h) | J^k. + \q } k.-\q) \ 

Q 1 J 1 G 1 + G 1 J 1 Q 1 - | 2E+{E2 _ _ E2+ _ + 2E _ {El _ E2 __ ^ j , (A.4) 



where we have retained the ie terms in the G\ propagators. As Eq. ( |A.3|) shows, these are 



not the correct ie factors, and they should be dropped immediately by taking the e — > limit 
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(as we are instructed to do). Dropping them gives a result identical to ( |A.3|) . Were we to 



(incorrectly) retain the ie's in Eq. ( A4) and expand the denominators into a principal value 
term and a delta function, the resulting delta function contributions to Eq. ( |A.4| ) would not 
cancel, giving an incorrect contribution to the current of the type Qi Jf Q\. We agree with 
Kvinikhidze and Blankleider that this contribution is spurious. It has not been included in 
any previous applications 0- ||, and is eliminated by taking the e — > limit [or simply 
dropping the ie terms from Eq. ( |A.4j )] after the contour integration has been carried out. 



APPENDIX B: GAUGE INVARIANCE FOR TRUNCATED CURRENTS. 

In this appendix we verify the gauge invariance of the truncated currents introduced in 
Section IV. 

First, for the purpose of further discussion it is convenient to split the divergence of 
the total untruncated current ( |3.6(J| ) into two parts corresponding to the divergences of the 
effective currents JfAofT an d Je Xe &, introduced in (|3.51|) and generated by the one-particle 
Ji A and the interaction Jg x Bethe-Salpeter currents, respectively. In particular 



QAeS = Qi hta), ] - [ ei (q), U}-U [eM + e 2 (q), A^] U Q 



'A ■ 



<l»JLe S = Qi [ei(q) + e 2 (q), U] + U [e x {q) + e 2 (q), A 9l ] U Q x 



(B.l) 
(B.2) 



The relation (B.l) follow from identities ( 3.38| - 3~42|) and its derivation can be repeated 
without any modification for the corresponding truncated effective currents. In deriving 
( B.2|) one has to use the quasipotential equation Q3.13j ) and more care is needed to get the 
divergence for the truncated J^ eff . 

Let us now consider the truncation by the rank N of the quasipotential U^ N \ as defined 
in eqs. ( |6.1| - |6.111 ). Using again the identities ( |3.38| - |3T42| ) one gets 



<b[ •>& + J&a + & ) = Gi(fc,(?), G- 2 1 } - [eM, V] + [eM + eM, V) ) Q x 

Q 1 ([e 2 (q),G^ + U^])Q 1 , (B.3) 



and repeating the derivation of ( |B.1| ) for truncated quasipotential with N > 1 

/ N ~ 1 \ 

Q*J$8& = Qi (-MQ), U {N) ] - E U {N - M) [eM + eM, A Pl ] ) Q 1 . 

V M=l J 

For the corresponding J^ls as gi ven by (|6.10| - |6.11|) we get 

Q„& = Qi(u^A 9l Ve(q) - U^A 9l e(q)V - e(q)VA 9l U^ 

N-2 v 

+Ve(q)A 9l U^ + U^ m ^ A gi [e(q), V]Ag 1 U^)Q 1 



(B.4) 
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<2i [e(g), U] + Ve(q) A 9l U^ + £ l/C^efe) A^C/W 

^ M=l 

N-2 v 

/ Ar ~ 1 \ 

Qi f[e 1 (g)+e 2 (g), [/] + E [/(^^(g) + e 2 (g), A^]C/W ) Q x , (B.5) 



M=l 



where we introduced the shorthand notation e(g) = ei(g) +e 2 (g) in intermediate steps. The 
derivation is valid for iV > 1, though for N = 2 some summations are empty. Clearly, the 
sum of O and flBJ) gives fl6T3l) . 

This derivation can be repeated for the case of two interactions defined by Eqs. ( |6.14| - 
|6.27|) . In this case one has to inspect the bounds of the summations carefully when the 
quasipotential equation is used, since the summations contain [/(°'°) = and therefore 
terms like ^Aftt/W^ 2 ' should be treated separately if M 1 = M 2 = 0. 

Finally, for the truncation by the number of exchanged mesons, as defined by Eqs. ( |6.2S - 
6.40| ), we obtained exactly as before 



n-l 



%J& = Qi l-lei(q), U^} - Y U {n - a) [ ei (g) + e 2 (g), A 9l ] 17« ) Q x 



a=l 



(B.6) 
(B.7) 



where (p3.7| ) is valid for n > and the sum does not contribute for n — 1. Both currents 
J e ( " } ^ from and flOOP can be considered at the same time and we get for n > the 

divergence 



n— 1 



a=l 



n-l 



_ Y U {n ~ a) A 9l V [a) 



a=l 



e(g) 



+ Y V (n - a) e(q)A gi U {a) - Y E U {n ~ a ~ h) Ag 1 V i - h) e{q)Ag l U { - a) 

a=l 6=1 a=l 

n— 1 n— 2 

- X] C/ (n ~ a) A#ie(g)y (a) + X] E ^ (n " a ~ 6) A^e(g)y (6) A^?7 (a) )Q! 

a=l 



=1 a=l 



n-l 



n-2 



+ £ V^e(q)A 9l U^ - Y 



0=1 

n-l 



a=l 



n— a— 1 



E C/( n -°- 6 )A^F^ 



6=1 



-(?)A^ (a) 



n— 2 n— 1 



EC/ (n-a) Aae(g)y (a) + i: E U^Ag 1 e{q)V^Ag 1 U^)Q 1 



a=l 6=1 c=6+l 

(n— 1 
Kg), f/ (n) ] + E f/ (n " a) e(g)A^ 1 f/ (a) 
a=l 

n— 1 n— 1 



c-1 



E U^ a) A 9l e(q)VW + E f/ (n ~ c) A^e(g) E V^A 9l U^ Q 1 



a=l 



c=2 



6=1 
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n-1 



Qi [e(q), UM] + £ £/ M [e(g), A^Jl/M & , (B.8) 



a=l 



where again some sums are empty for n — 1, 2. The sum of ( |B.7| ) and ( [B.8| ) yields 
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